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Abstract 



Molecular dynamics computer simulations which employ the embedded- 
atom potential show that nanowires of gold exist as multishelled structures. 
We simulate double-walled gold nanowires and calculate the capacitance of a 
finite nanometer-size cylindrical capacitor. For the sizes for which multishelled 
nanowires appear in simulations we find the capacitances below one attofarad. 
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Finite-size systems of atomic scale dimensions have attracted great attention. Interest 
is focused on fundamental aspects, as well as on applications. The remarkable physical 
properties of nanostructures allow one to design nano-scale electrical and mechanical devices. 
The continuing miniaturization of engineering devices leads to a technological revolution. 
Research on the properties of cylindrical nanostructures of various materials, ie., nanowires, 
has been especially active field. 

Multishelled nanostructures were produced for carbon clusters and wires [0,0, as well 
as for wires of WS 2 , MoS 2 , and NiCl 2 @. Gold nanowires are fabricated by Scanning 
Tunneling Microscopy (STM) and electron-beam litography |J. In recent Molecular 
Dynamics (MD) computer simulation multishelled cylindrical nanostructures of gold were 
obtained @. Wires with radii around a nanometer and a length/diameter ratio between 
1 and 3 were studied. It was found that multiwalled structures exist for longer of these 
nanowires, i.e., for a length/diameter ratio of 2 and 3. The formation of shells was also 
found in a jellium model calculation for finite sodium nanowires ||. This was confirmed 
by the conductance measurements ||. We simulate here double- walled gold nanowires and 
study the capacitance for this geometry. 

It is well known that an application of a many-body potential is necessary for accurate 
description of metallic bonding by the MD simulation method [TIJ . We used here (as in Ref . 
0) an embedded-atom potential for gold which produces a good agreement with available 
experimental results for bulk, surfaces, and nanoparticles fLl |. A time step of 7.14 x 10~ 15 



s was employed in simulation. The temperature was controlled by rescaling the particle 
velocities. We started from an ideal fee (111) structure at T = K, and included in the 
cylindrical MD box all particles whose distance from the nanowire axis was smaller than 1.2 
nm. The length of the sample of 689 atoms was 6 layers. We also studied a nanowire of 
9 layers and 1032 particles. The samples were first relaxed, then annealed and quenched. 
To prevent melting and collapse into a drop, instead of usual heating to ~ 1000 K used in 
MD simulation of gold nanostructures our finite nanowires were heated only to 600 K. Such 
a procedure, as in Ref. 0, gives the atoms a possibility to find local minima and models 
a constrained dynamical evolution present in fabricated nanowires. The structures were 
analyzed after a long MD run at T = 300 K. Figure 1 shows the shape of the MD box for 
a nanowire of 689 atoms after 3.55 ns of simulation. A top view of the whole cylindrical 
MD box for this nanowire is shown in Fig. |2|. While the presence of a double-walled and 
internally filled structure is obvious, the walls are still not completely homogeneous. Several 
atoms remain about the walls. Simulation shows that an almost identical structure exists 
up to 10 6 time steps, i.e., 7.1 ns fl2H. We also found a double- walled structure for nanowire 



of 1032 atoms. These structures suggest an application of gold nanowires as cylindrical 
capacitors. Therefore, we calculated the capacitance of a nanometer-scale capacitor for the 
geometry shown in the inset of Fig. |[ 

The calculation starts from the electrical potential 0(r, z) in the space between two 



coaxial cylinders of radii r\ and r 2 (r x < r 2 ), and finite length L in the z direction |L3 



where v k = kir/L, r\ <r < r 2 , and < z < L. The parameter V is the potential difference 
applied between two cylindrical plates, whereas h{x) and Ki(x) (here and below) are the 
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modified Bessel functions of the order i (T^|. The capacitance C = (ri, r 2 , L) of the structure 
is obtained from 

\ / r=r± 

where eo = 8, 854 x 10~ 12 F m~ l is the electric permittivity of vacuum. The result is 

C = 16r e V 1 ^0(^1)^1(^1) + Ji(^ri)ir (z/ fc r 1 ) 
1 ° fc= i~^ 5v k Io(v k ri)Ko(v k r 2 ) - h[v k r 2 )K^{v k r x ) ' 

Using equation (3) we calculated the capacitance for several sizes typical for structures 
obtained in MD simulations. The results are presented in Fig. |3] and Table 1. 

In summary, molecular dynamics computer simulation based on the well-established po- 
tential shows that gold wires of nanometer dimensions form double-walled structures. Mo- 
tivated by this result we calculated the capacitance of a finite nanometer-scale cylindrical 
capacitor. Our calculations of the capacitance is based on the classical electrodynamics. 
However, on very small scales the quantum corrections appear. This is the point we are 



now pursuing |15[]. It is known that typical inhomogeneous electric fields used in the po- 
larizability measurements for metallic clusters are on the order of 10 5 Vcm~ x and field 
gradients are ~ 10 5 Vcm~ 2 JT^|. These fields induce a dipole moment and slightly deform 
clusters. Therefore, to strictly preserve the shape of a nanometer-scale capacitor weaker and 
homogeneous fields should be applied in an electronic device. An STM could be used to 
manipulate gold atoms and assemble small nanowires j|,[5|] . The STM experiments show that 
gold nanostructures can sustain large fields of up to 1 V 0. Wires produced by electron- 
beam litography will approach the smallest nanometer-scales discussed here. It is possible 
to fill the space between two cylindrical electrodes with various insulators and to obtain 
nanometer-scale capacitors with different capacitances. Similar filling of carbon nanotubes 
with metals was recently realized [|1^,|T^]. Moreover, fabrication of metallic nanowires us- 
ing carbon multishelled nanotubes as templates could be used to produce nanometer-scale 
cylindrical capacitors not only of gold, but also of other metals. 
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TABLES 

TABLE I. Capacitance C of a nanometer-scale capacitor with an external radius T2 = 1.2 nm, 
an internal radius r\ = 0.9 nm, for several lengths L. 



L(nm) 


2.0 


2.2 


2.4 


2.6 


2.8 


3.0 


C{aF) 


0.336 


0.374 


0.413 


0.452 


0.490 


0.529 
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FIGURES 



FIG. 1. Atomic positions (side view) for a nanowire with length 2.6 nm and an external radius 
of 1.2 nm. 

FIG. 2. MD trajectory plot of a top view for a gold nanowire shown in Fig. 1. The whole 
thickness of the wire along its axis is shown after 3.55 ns of simulation. This nanowire consists of 
the two coaxial near walls and a large filled core. 

FIG. 3. Capacitance of a capacitor as a function of an internal radius n, for L = 2.6 nm, and 
an external radius T2 = 1.2 nm. Inset: the geometry of a nanometer-scale cylindrical capacitor. 
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This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/cond-mat/0005276vl 



